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Neutrino masses from operator mixing 
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Abstract 

We show that in theories that reduce, at the Fermi scale, to an extension of 
the standard model with two doublets, there can be additional dimension five 
operators giving rise to neutrino masses. In particular there exists a singlet 
operator which can not generate neutrino masses at tree level but generates 
them through operator mixing. Under the assumption that only this operator 



appears at tree level we calculate the neutrino mass matrix. It has the Zee 



, mass matrix structure and leads naturally to bimaximal mixing. However, the 



maximal mixing prediction for solar neutrinos is very sharp even when higher 
order corrections are considered. To allow for deviations from maximal mixing 
O • a fine tuning is needed in the neutrino mass matrix parameters. However, 

this fine tuning relates the departure from maximal mixing in solar neutrino 
O . oscillations with the neutrinoless double beta decay rate. 



PACS numbers: 14.60.Pq, 14.60.St, 14.80. Cp, 12.60.Fr 



The simplest model for neutrino masses is based on the seesaw mechanism In the 

seesaw mechanism the standard model (SM) is enlarged with singlet right-handed neutrinos. 
Then, a Dirac mass term, Md, mixing left-handed and right-handed neutrinos is possible. 
In addition, since right-handed neutrinos do not carry any gauge charge, they can have 
a Majorana mass, Mr, without compromising the gauge symmetry. If the right-handed 
Majorana mass term is very large, as expected for a singlet mass term, very light Majorana 
neutrino masses for left-handed neutrinos are obtained through the diagonalization of the 
full mass matrix of neutral fermions, m u = M^/Mr, thus justifying the small size of neutrino 
masses. Since the Dirac mass term is proportional to the standard Higgs vacuum expectation 
value, this mechanism provides masses which £1X6 771 ]y ~ A|i/A, being A^ the Fermi scale and 
A the lepton number breaking scale. This type of behaviour is much more general and, in 
fact, many neutrino mass models can be cast into this form. This can be understood in the 
following way: if the SM is just the low energy effective manifestation of some underlying 
theory, the effects of new physics can be represented by a series of gauge invariant operators 
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containing the SM fields with higher dimension operators suppressed by powers of the scale 
of new physics [f§-f|. At low energies, the most relevant operators will be those with the 
lowest dimension, that is dimension five operators. One can easily see that there is only 
one gauge invariant operator of dimension five one can build with the field content of the 
standard model 

^seesaw = ~j{£Fr£)(^T(f) , (1) 

where £ is the standard left-handed doublet of leptons, £ = iT 2 £ c , l c = C£ ( C is the charge 
conjugation operator), ip is the Higgs doublet and (p = iT 2 (p*, f are the Pauli matrices in 
SU(2) space, F is a complex symmetric matrix in flavour space (577(2) and flavour indices 
have been suppressed) and A is a scale related to the scale of new physics. It is clear that 
this Lagrangian does not conserve generational lepton numbers, but in addition it does not 
conserve the total lepton number, which is violated in two units. In the SM, lepton number 
is conserved because the requirement of renormalizability and the small particle content (no 
right-handed neutrinos, no triplet scalars, etc) but as long as the spectrum is enlarged there 
is no strong reason for lepton number conservation. Operator (|l|) makes this statement 
explicit. Therefore, this operator will be generated in any extension of the SM that does 
not conserve lepton number. 

When the Higgs develops a vacuum expectation value (VEV), operator ([]]) will give rise 
to a neutrino Majorana mass matrix given by 

M„ = Fj, (2) 

with v = ^v 9 ^) — 174 GeV, the SM Higgs VEV. If we take the largest eigenvalue of F 
to be of order 1 and use the laboratory bound on the r-neutrino mass, m„ T < 18 MeV, 
we find that A > 10 6 GeV. Should one take the value suggested by atmospheric neutrino 
data, m Vr pa 0.06 eV, one would obtain A pa 5 x 10 14 GeV, a scale which is close to the 
unification scale: physics at very high energy scales could affect low energy physics in a 
measurable way through neutrino masses. However the relationship between A and masses 
of new particles can be quite different from the naive expectations. For instance, it seems 
that the Lagrangian of eq. (P is generated by the exchange, among the leptons and the 
Higgses, of a scalar triplet, Xi with hypercharge 1. Then, 1/A ~ fi/mi with /x the trilinear 
coupling of the triplet with the Higgs doublets. However, this is not the only possibility. In 
fact, eq. ([]]) can be identically rewritten, after a SU(2) Fierz transformation, as 

^seesaw = -"(F^)F(^ , (3) 

which suggests the exchange of a neutral heavy Majorana fermion; then A should be the 
mass of that fermion. Indeed, the seesaw mechanism described above implements naturally 
this possibility. 

Given the full generality of this mechanism, which naturally relates the smallness of 
neutrino masses to the new physics scale, it is quite reasonable to try to fit the different 
neutrino mass models into this description. However, this is not always possible: on one side 
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it can happen that neutrino masses are generated only by operators with higher dimensions 
(for a recent analysis of the different possibilities see 0) on the other side models in which 
neutrino masses are generated through radiative corrections are also difficult to fit in the 
simplest scheme. 

In this paper, we want to show that the uniqueness of the effective seesaw mechanism 
can be relaxed a bit if new fields are allowed at the electroweak scale, in particular in models 
with two light doublets. This opens the door for a new class of effective seesaw mechanism 
in which the light neutrino masses are generated radiatively. This fact will allow us to 
lower the lepton number breaking scale by several orders of magnitude. In addition, as 
we will see, this mechanism predicts a very particular form for the neutrino mass matrix 
which seems well suited for explaining both atmospheric and solar neutrino data. This is 
not strange since a particular realization of this mechanism is the Zee model ||,^| and its 
variations. For instance, models with spontaneous symmetry breaking of lepton number by a 
doublet |nj . Another interesting variation are models in which spontaneous lepton number 
breaking is induced by a hyperchargeless triplet fTl,l!^|. They have the interesting property 



that although the triplet feels gauge interactions it does not contribute to the invisible Z 
decay width and that the lepton number breaking VEV could be at the electroweak scale (no 
bounds from red giant cooling by Majoron emission). Models with hyperchargeless triplets 
have also been discussed recently |13],14] in connection with the results of last standard model 



fits. Variations solving the strong CP problem can be found in |15|,[16|. All of these models 
predict a Zee type neutrino mass matrix (for a comprehensive review of extensions of the 
Higgs sector of the SM see |17|,[18|). Recent fits to neutrino data |19H22|1 suggest some type of 



bimaximal mixing which can be accommodated naturally in this type of models p3|-|26|j . This 



observation has boosted again the interest of models with Zee type neutrino matrix [27-33 



In the same way that a variety of see-saw mechanism models can be described as a single 
effective operator it would be interesting to see if this class of models and perhaps other type 
of models with two light doublets can be described at low energies with just a few operators. 

We will assume that the low energy (Fermi scale) theory is just the SM model, with 
no right-handed neutrinos, supplemented by an additional doublet. We will denote the two 
doublets as <pi and (p 2 . Then, in principle, nothing forbids that both doublets couple to the 
two types of quarks and to the leptons. However, it is well known that this, in general, will 
lead to neutral current flavour changing problems pi |. Therefore, for the moment we will 
consider that only one of the doublets, tpx, couples to the fermions. This can be archived 
naturally by assigning an additional conserved charge to the doublet tp2, lepton number for 
example. Of course, this additional charge should be explicitly broken in the Higgs potential 
in order to avoid the appearance of a Goldstone boson once the doublet <p 2 acquires a VEV 
(models in which lepton number is broken spontaneously by a doublet [|l(],|35],[3(J are excluded 
by the invisible Z width measured at LEP). Therefore, for the moment we will assume the 
SM Yukawa couplings to leptons 



Cy = £YeRtpi + h.c. 



(4) 



Again, SU(2) and flavour indices have been suppressed and Y is a 3 x 3 complex matrix 
in flavour space. However, because no right-handed neutrinos have been introduced in the 
model and because only (fi couples to leptons, it can be chosen, without loss of generality, 
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as diagonal with all its elements real and positive. As in the SM, neutrinos remain massless 
because there are no right-handed neutrinos and because lepton number is conserved. 

Now we will assume that this model is just the low energy manifestation of a more 
complete theory which will only shows up at higher energies. As discussed above, if the 
scale of new physics is high enough its effects can be be parametrized by operators with 
higher dimensionality: 

C e ff = Ahsm + j^Cx + ^2^2 + ■ ■ • • (5) 

Here ^2Hsm represents the renormalizable Lagrangian we just outlined which is a minimal 
extension of the SM containing two doublets. 

If the theory at the Fermi scale contains two doublets, one can write four independent 
dimension five operators^ 

C ^ = -lUmre^r^), £ T2 = - l -\(J.T 2 Tt)^T V2 ) , (6) 



C T = -—(ITTiW^) , (7) 

and 

c 3 = -~(Jse)(^i Pl )- (8) 

Operators Ct x , Ct 2 can be excluded by the same symmetry used to forbid Yukawa cou- 
plings of the doublet ip 2 to the fermions. For instance, one can assign lepton number to (p 2 in 
such a way that it does not have Yukawa coupling to fermions while allowing the couplings 
Ct and Cs- L({p 2 ) = —2 will do the job and will forbid operators C^, Ct 2 Therefore, in the 
following we will only consider operators Ct and Cs- 

It is important to notice that the operator in (H) does not exist with only one doublet 
because the singlet coupling of two scalar doublets is antisymmetric (in SU (2) components 
it is just eijif 2 iifij). In addition, one can see that, due to Fermi statistics, the 3x3 matrix 
in flavour space, T, is symmetric while the singlet coupling S is a complex antisymmetric 
matrix. 

When the Higgs fields develop a charge conserving VEV, operator Ct gives rise to a 
neutrino Majorana mass. However, the singlet operator Cs does not give any mass to the 
neutrinos because the product of the two doublets has hypercharge 1, and a singlet with 
hypercharge 1 has also charge 1 (at difference with the triplet combination (|7|) which has 
neutral components). For this reason the singlet operator does not seem very interesting at 
first sight. However, a very interesting situation arises when, for some reason (limited particle 



x In general one also expects higher dimension operators which could contribute to interesting 
processes as [i-e conversion in nuclei, [i — > &y, etc. For instance, when a charged scalar is integrated 
out all those operators appear at one loop |}7). 
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content of the full theory), only operator @ arises at tree levelQ. As commented before, 
it cannot give rise to neutrino masses after spontaneous symmetry breaking. However, one 
expects that renormalization effects will mix all operators with the same dimensionality and 
the same quantum numbers. Therefore, (§) and (]?]) will mix under the renormalization 
group and operator ([!]) will be generated at one loop even if it did not appear at tree level. 
In fact, one can easily see, by computing the diagram in fig. [I] (and the crossed diagram) 
and taking the divergent part, that the matrix T obeys the following renormalization group 
equation (RGE) 

Here, the dots represent extra contributions to the renormalization group of the T matrix 
which are proportional to the T matrix itself and, therefore, cannot generate it if it did not 
exist at some scale. This RGE is peculiar with respect the RGE we are used to see for the 
Yukawa couplings, in both the SM and the MSSM, in what there is a piece that does not 
depend on the T matrix. In both the SM and the MSSM there are several chiral symmetries 
broken only by the Yukawa couplings that ensure that the RGE of those couplings should 
transform in a covariant way with respect to those symmetries. This ensures, in this type 
of theories, that fermion masses cannot be generated through radiative corrections. Equa- 
tion @ is not of this type and, therefore, even if the coupling T did not exist at the scale A 
at which the operators were generated, it will arise through operator mixing. It is very easy 
to integrate eq. (|9]) by keeping only the leading logarithm (a more sophisticated integration 
can be performed by taking into account the running of the standard Yukawa couplings, the 
running of the S coupling itself and the extra couplings we neglected in eq. (|), however, 
this effect is higher order in the couplings and since the couplings are small we expect a 
small effect). The result is that 

T(m z ) » JL (SYY^ + Y*Y T S T ) log (^) + T(A) , (10) 

where we have identified mz with the Fermi scale. We will assume that T(A) is not generated 
at tree level. Of course T(A) could also pick up contributions at one loop (or higher loops). 
To compute them one would need to know the details of the full theory in which our effective 
theory is embedded. However if A is large enough T(mz) will be dominated by the model 
independent logarithmic piece in eq. ( |10D which can be computed in the effective theory. So, 
as a first approximation we will assume T(A) ~ and further on we will keep T(A) only 
when the logarithmic pieces vanish. 

After spontaneous symmetry breaking, if both (pi and y?2 develop a VEV, operator ([?]) 
will give rise to a neutrino mass matrix for the lefthanded neutrinos given by 



2 An example of this situation is provided by the Zee model in which a charged singlet, h + , and an 
extra doublet are introduced with couplings £ f£h + and \ihr(fyp\, then, for a heavy h + one finds, 
at tree level, that S/ A = /mi while the triplet operator cannot be obtained at tree level. 
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M " " Ual3 Wf ( Si ' yt + r " yTST ) '° g (x) A ' (U) 

where V\ = (^(pf^, v 2 — (v 9 ^) are the VEV's of the two doublets and tan/? = v 2 /vx. 
The seesaw structure is apparent in the last term. The other factors, however, are also 
important. First, the neutrino mass matrix comes naturally proportional to the mass of the 
leptons squared which gives an important suppression since lepton Yukawa couplings are 
small. Second, it contains the standard loop suppression factor 1/(Att) 2 and, third, the ratio 
of VEV's of the second doublet v 2 to the standard doublet VEV, tan/5 = v 2 /vi, can give an 
additional suppression factor. All together, one can achieve with this mechanism the same 
neutrino masses one would obtain in a standard seesaw mechanism with a A which is at 
least 6 orders of magnitude smaller. This could put the scale of the new physics responsible 
for neutrino masses at the reach of the next generation of accelerators if the ratio of VEV's 
v 2 /vi and/or the largest Sij are very small. However, perhaps the most interesting aspect of 
eq. ( JTT| ) is the structure of the mass matrix, inherited from the antisymmetric structure of 
the singlet coupling S: if we choose for the Yukawa coupling of the leptons a diagonal form 
we find 



M v « — m 
with x e = ml/m 2 , x^ = m 2 Jm 2 and 



S e/1 (x e - x^) - 1) | (12) 

\ S eT {x e - 1) Sprfap - 1) 



rn 
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m = tan/3 r log . 

A(47r) z \mz / 

The structure of this mass matrix is very interesting since all diagonal elements are zero and 
contains only three free parameters (the three elements of the antisymmetric matrix S). It 
is convenient to rewrite it as 




M v ~ m efl 77V ■ ( 13 ) 
\ m eT m^r / 

This form of mass matrix has been considered recently in order to fit both atmospheric 



and solar neutrino data |23l- f2q| . Let us review some of the results: in order to explain 
solar neutrino data in terms of oscillations one needs a mass squared difference, A s which is 
10~ n eV 2 < A s < 10~ 5 eV 2 . A global analysis including the results of SNO suggests mixings 
close to maximal. On the other hand, in order to explain atmospheric neutrino data one 
needs a mass squared difference, A a w 3 x lCT 3 eV 2 and also a very large mixing. The 
particular structure of the obtained mass matrix (it is traceless) implies that the sum of 
the eigenvalues is zero a fact which constraints the possible solutions. An analysis of the 
different possibilities in terms of this mass matrix has been carried out in [f23| -|26| where it 
has been shown that only the case with m e[1 ~ m eT and m eil , m er is acceptable. This 

naturally predicts maximal mixing for solar oscillations which, after SNO, seems to be the 
only viable possibility. In fact the Zee mass matrix predicts, in this case P7j, 
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sin 2 26, 



16 \A n 



(14) 



This is a very strong prediction which is compatible with the LOW and the VAC solutions of 
the solar neutrino problems. The LMA solution, which right now seems to be the preferred 
solution, would be marginally compatible with this prediction. However, it is important to 
notice that in general we expect corrections to the Zee mass formula. By requiring that 
only one doublet couples to fermions we have taken the most restrictive (and predictive) 
possibility. By relaxing this assumption one can perfectly fit also the LMA region in solar 
neutrino parameters |j28 |. But even in the restrictive case in which only one doublet couples 



to fermions one can also take into account subdominant contributions. In fact, on general 
grounds we expect modifications to the particular form of the mass matrix with all the 
elements in the diagonal vanishing. There is no symmetry that enforces this structure, 
therefore one expects that at some point this structure will receive corrections. This happens, 
for instance in the Zee model where diagonal entries are generated at two loops []3C|| . For 
our purposes we can include all those corrections in the initial contributions at the scale A, 
that is, by considering a non- vanishing T(A) which we could take as a general symmetric 
matrix. Given the bimaximal mixing required by the data it is natural to parametrize the 
neutrino mass matrix as follows 



M v « 



m 



73 
\T2 



^ Tee 1e/i 7er 
1e[i 7mm 7/-4T 
\ 7er 'Jiit ^tt 



(15) 



where, without loss of generality we can take 7 er = — 7 eM and choose one of the 7's equal to 
one defining in this way the normalization of e. We will assume that the first term, generated 
through running and containing the logarithmic enhancement, is the dominant one and gives 
the scale of atmospheric neutrinos. The term proportional to e is subdominant and cannot be 
computed in the effective theory. Its explicit form depends on the details of the underlying 
theory and it could contain one loop contributions not enhanced by large logarithms or/and 
higher loop contributions, depending on the underlying theory. For simplicity, we will also 
take all the 7's as real. Notice that the mechanism suggested in this paper only gives a leading 
order contribution with zeros in the diagonal. The extra structure assumed in the leading 
mass term ( m efl m m eT and m^ T « 0) must be imposed with some additional symmetries 
which, however, are not unnatural in this type of models ||. The modifications introduced by 
the term proportional to e can be qualitatively very important. One of the most interesting 
predictions of the leading order mass matrix is that there is no neutrinoless double beta 
decay (NDBD) : since the NDBD amplitude is proportional to (M u ) ee it is obvious that 
a mass matrix with zeroes in the diagonal forbids NDBD. However any correction to the 
leading order mass matrix introducing diagonal entries will lead to NDBD. On the other 
hand, corrections to the leading order mass matrix are necessary in order to accommodate 
solar neutrino mass differences and small departures from maximal mixing in solar and/or 
atmospheric neutrinos. It is important to check how these corrections could modify the 
sharp predictions of the model. 

By diagonalizing the mass matrix at second order in e one easily obtains that 
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A a w ml , (16) 
sin 2 2# a « 1 - - (l6 7eM + (7 W - Irrf) e 2 , (17) 
A s « m 2 (2 7MT + 2 7ee + 7ftU + 7rr )e , (18) 

sin 2 2fl s « 1 - -L + i ( 7ee (2 7MT + 7w + 7rr ) - ( 7wi - 7rr ) 2 ) e 2 . (19) 

Since A s /A a is known to be small, the natural prediction of the model is just maximal 
mixing for the two angles to a very good degree of precision. Departures from maximal 
mixing are allowed naturally for the atmospheric mixing angle since A s does not depend 
on 7eA1 which controls the atmospheric mixing (for 7w , = 7rr ), therefore it can be made 
large without any conflict with solar neutrino data. However, to accommodate deviations 
from maximal mixing in the solar neutrino parameters is more delicate. One needs to make 
A s small while keeping the last term in eq. (fT9"D relatively large. This is clearly unnatural. 
For instance, in the case in which 7w< = 7rT = one obtains A s /A a w 2 ( 7mt + 7ee ) while 
sin 2 29 s ~ 1 — (A s /A a ) 2 /16+ 7Mr7ee . Therefore, to obtain a very small A s /A a while keeping a 
sizeable contribution to sin 2 29 s one would need to fine tune the couplings in such a way that 
7ee « — 7mt . Although it is not natural, this possibility might be interesting because perhaps 
it is the only way to accommodate LMA within this scheme and because if it is realized in 
nature it will link the deviation from maximal mixing in solar neutrino oscillations with 
neutrinoless double beta decay: as commented before, the NDBD amplitude is proportional 
to the so-called effective neutrino mass (m u ) = U^m^ = (M v ) ee . 

In this work we have shown that in models that contain in its low-energy (Fermi scale) 
effective theory two Higgs doublets there are four independent dimension five gauge invariant 
operators which violate lepton number. Three of them, which couple leptons to doublets 
in the triplet channel, generate masses at tree level when the doublets acquire VEV's. The 
other operator, which couples leptons and doublets in the singlet channel, cannot generate 
masses at three level. However, loop corrections mix all operators under the renormalization 
group and, therefore, the singlet operator also gives rise to neutrino masses when the doublets 
acquire VEV's. 

Under the assumption that only the singlet operator is generated at tree level, at the 
scale of new physics, one can compute the induced neutrino masses at one loop. Neutrino 
masses are suppressed by several factors, the loop factor, the masses of charged leptons and 
the ratio of the two doublets VEV's, therefore allowing for a new physics scale several orders 
of magnitude lower than the one needed in tree-level mechanisms for neutrino masses. 

Because the structure of the singlet effective operator, which necessitates antisymmetric 
couplings in flavour, the obtained mass has the structure of the Zee mass matrix with zero 
entries in the diagonal. This structure naturally accommodates bimaximal mixing and is well 
suited for explaining both solar and atmospheric neutrino data. However, in the simplest 
scheme it is very difficult to accommodate the LMA solution for solar neutrinos because it 
gives a sharp prediction for sin 2 29 s = 1 once one takes into account the mass differences 
needed for solar neutrinos. 

We have investigated the possibility to accommodate sin 2 29 s ^ 1 in this scheme by 
considering non-leading contributions to the mass matrix since, on general grounds, one 
expects to generate non-zero entries for the diagonal elements of the mass matrix. We found 
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that the prediction sin 2 26 s = 1 is quite stable and, only by fine tuning the parameters of the 
mass matrix, it is possible to accommodate at the same time s\vl26 s 7^ 1 and A s /A a <C 1. 
This fine tuning, even though unnatural, has the interesting property that relates the rate 
of neutrinoless double beta decay to the departure from maximal mixing in solar neutrino 
oscillations. 

This work has been funded by CICYT under the Grant AEN-99-0692, by DGESIC under 
the Grant PB97-1261 and by the DGEUI of the Generalitat Valenciana under the Grant 
GV98-01-80. 
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FIGURES 




FIG. 1. Diagram contributing to the mixing of the singlet and triplet operators 
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